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Abstract 

The effects of an externally applied one-dimensional periodic potential on the freezing/melting 
behaviour of two-dimensional systems of colloidal particles with a short-range attractive interaction 
are studied using Monte Carlo simulations. In such systems, incommensuration results when the 
periodicity of the external potential does not match the length-scale at which the minimum of 
the attractive potential occurs. To study the effects of this incommensuration, we consider two 
different models for the system. Our simulations for both these models show the phenomenon of 
reentrant freezing as the strength of the periodic potential is varied. Our simulations also show 
that different exotic phases can form when the strength of the periodic potential is high, depending 
on the length-scale at which the minimum of the attractive pair-potential occurs. 

PACS numbers: 82.70.Dd, 64.70.Dv 
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I. INTRODUCTION 

In the pioneering experiments of Chowdhury, Ackerson and ClarkjlJ on laser induced 
freezing, a two-dimensional monolayer of colloidal particles in the liquid state was subjected 
to a laser intensity pattern periodically modulated along one direction. They found that if 
the wave-vector of the modulation is tuned to the wave- vector at which the liquid structure 
factor peaks, a triangular lattice with full two-dimensional symmetry results for laser intensi- 
ties above a threshold value. In a later experiment Wei et al P] observed that this triangular 
lattice melts if the strength of the laser field is increased further. Although such a 'reentrant 
melting' phenomenon was observed in Monte Carlo studies [^] of colloidal particles inter- 
acting via Derjaguin-Landau-Verwey-Overbook(DLVO) potential, later simulations of the 
same system by Das et al [jj did not show the re-entrant phase. However, by extending 
the Kosterlitz-Thouless-Halperin-Nelson-Youngi6] theory of defect-mediated melting in two 
dimensions to the case where an external periodic potential is present, Frey et a/jjj argued 
that for short-range interactions there will indeed be a re-entrant transition to a liquid state. 
Later experiments (^J and simulations [l| of charged colloids seem to confirm the occurrence 
of re-entrant melting. Meanwhile, numerical studies 
periodic potential on a system of hard disks have also shown the occurrence of re-entrant 
melting. 

All the above results correspond to the case when the commensurability ratio p = \^3a/2d 
(where a is the mean particle distance and d is the period of the external periodic poten- 
tial) has the value of 1, i.e every potential trough is occupied by particles. In a recent 
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10, hj of the effect of the external 



experiment 



121 ] corresponding to p = 2, a new phase - the 'locked smectic state' - was ob- 



served, with the crystalline state being found to melt via this new intermediate phase. This 
observation is also in qualitative agreement with theoretical predictions j^]. 

No study has been reported as yet for the case where there is a possibility of incommen- 
suration between the periodicity of the external laser field and some length-scale inherent to 
the two-dimensional colloidal system. In the present work we study just such a case. Specif- 
ically, we look at the effect of the external laser field on a monolayer of colloidal particles, 
which interact via a short-range attraction apart from the usual hard-core repulsive interac- 
tion. Such a short-range attraction is known to arise when hard sphere colloidal particles are 
mixed with polymers or smaller colloidal particles, giving rise to an effective attraction ["K^ 
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between the larger particles, called depletion interaction. 



context of colloidal particles confined between two walls 



t has also been suggested in the 



15] 



14, 13, where it is claimed that an 



attractive minimum in the effective interaction between a pair of colloidal particles arises 



when they are close to the walls 



161 ] . Effects due to incommensuration should occur if the 



length-scale at which the attractive part of the interaction has its minimum is different from 
the periodicity of the externally applied laser field. Novel phases can then appear in the 
system depending upon the relative strengths of the two competing potentials with different 
length scales. 

Recent!, GoUe « al Q studied the effect of an externa, potent., witn one-d—al 
periodic modulation on the liquid-vapour transition line of a three-dimensional mixture of 
hard-sphere colloidal particles and polymers. In their density functional calculations, they 
obtained a new stacked fluid phase, which consists of a periodic succession of liquid and 
vapour slabs. The presence of this new phase is a manifestation of the one-dimensional 
nature of the externally applied modulated potential in a three-dimensional fluid. Their 
simulations and calculations also showed that the density profiles exhibit a non-monotonic 
crossover when the wavelength of the modulation matches the hard sphere diameter. 

In the present work, we analyze the effects of incommensuration on the phenomenon of 
freezing/melting of colloidal particles in two dimensions in the presence of a tunable "sub- 
strate" or external potential with one-dimensional periodic modulation. For this purpose, 
we consider two different models that are described in section 2, and use Monte Carlo simu- 
lations to study their behaviour. In our simulations of these two models, we clearly see that 
for both low and high values of the external potential the colloidal particles form crystalline 
structures, whereas for intermediate strengths of the external potential the system is in a 
modulated liquid phase. We term this phenomenon re-entrant freezing. It is also observed 
that at high values of the external potential, different kinds of phases can occur depending 
on the length-scale at which the attractive part of the inter-particle potential has its mini- 
mum. Our simulation results are described in detail in section 3. The main conclusions of 
this study are summarized in section 4. 
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II. THE MODELS 



A. Hard-core particles with short-range attraction 

The generation of an effective attractive interaction between large hard-sphere colloid 
particles when they are mixed in a binary mixture with small hard-sphere colloid particles 
has been known for quite some time. In a recent study Castafieda-Priego et al [3| calculated 
this depletion interaction in a strictly two-dimensional binary mixture of hard disks. Using 
the expression of the depletion interaction derived in their work, we have constructed a model 
for a two-dimensional system of colloidal particles with short-range attractive interaction. 
The pair potential Ufrij) between particles i and j with centre-to-centre distance is given 
by: 
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Here, a is the particle diameter and U and rj are two parameters which can be used 
to tune the depth and width of the short-range attractive part of U(rij). This interaction 
potential has a minimum at = a, i.e. when the two particles touch each other. The 
strength of the interaction potential at this minimum is given by 
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In addition, we assume that a particle with co-ordinates (x, y) experiences an external 
periodic potential of the form 



27T 

V[x,y) = V e cos(—y) 



(2.3) 



In the above equation, the constant d is chosen such that for a number density of 
= N/(L x L y ) (N is the number of particles and L x , L y are the lengths of the sides of 
a rectangular sample), the modulation is commensurate with a triangular lattice with near- 
est neighbour distance a s = 1/(^0)5, i.e. d = a s a/3/2. Thus, for the binary mixture we 
are trying to model, the wave-vector of the external potential is commensurate with the 
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smallest reciprocal lattice vector (q = 27r/(a s y / 3/2)) of the triangular lattice that the large 
disks would form at that density in the absence of the smaller disks. Also, the form of the 
external potential is such that its troughs run parallel to the 

For the system under study, the important parameters are U = f3U min , j3V e , p = (pa 2 and 
7], where j3 = l/(fcsT), fee being the Boltzmann constant and T the temperature. While 
the attractive part of U{rij) would like to have the particles touching each other (with 
inter-particle separation a), the externally applied potential would like to have a density 
modulation in the y direction with periodicity d, resulting in the incommensuration. 



B. Soft-core particles with short-range attraction 

The second model we consider is one where the attractive potential is considered to be 
a Gaussian well, with the position of its minimum being incommensurate with the inter- 
particle separation. Thus, the pair potential U(rij) between particles i and j with distance 
Tij is given by 

U(r „) = l[(^(^) 2 5^2^_ J 4 exp{ _ B(r , j _ A)2} ], (2 . 4) 

where R is the radius of the colloidal particles with surface charge Z*e and k is the inverse of 
the Debye screening length. Similarly to the first model, in this case also a colloidal particle 
with co-ordinates (x, y) is assumed to experience an external periodic potential of the form 

2tt 

V(x,y) = V e cos(—y) (2.5) 
a 

As in the first model, the modulation of the external potential, d is chosen such that for 
colloidal particles with density = N/(L x L y ), the modulation is commensurate with a 
triangular lattice with nearest neighbour distance a s = 1/(^0)^, i.e. d = a s v^3/2. The 
parameter A in Eqn l2.41 which determines the position of the minimum of the attractive 
part of the potential, is assumed to be incommensurate with the triangular lattice. The 
parameter A determines the depth of the attractive well and B determines its width. 

This model was motivated in part by the suggestion [l^. Q] that in a two-dimensional 
monolayer of charged colloidal particles in an aqueous solution confined between two walls, 
with the 'bare' interaction between pairs of particles being given by the Derjaguin-Landau- 



Verwey-Overbook(DLVO) potential 
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an attractive well with a depth of the order of ksT 
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develops in the effective potential between a pair when the particles are near a charged wall. 
The depth and the position of the minimum are supposed to be strongly dependent on the 
distance from the wall, the effective charge on the colloid particles, the counter-ion density 
and the surface charge on the wall. However, this claim has been contested^]. In any 
case, leaving aside the question of a physical realization of the above model, its use tests the 
robustness of our results with respect to the details of the model interaction. 



III. SIMULATION AL DETAILS AND RESULTS 



A. Model 1 



We have carried out Monte Carlo simulations of a system of N particles of diameter 
a, interacting via the potential Ufcj) defined in Eqn l2.ll The particles are contained in 
a rectangular box of dimension L x x L y , where (L y /L x = v3/2), with periodic boundary 
condition being used for doing the Monte Carlo simulations. Most results reported here are 
for simulations done for N = 1600 particles. In order to check finite size effects, we will also 
discuss results for iV = 1024 and 900. 

For a system of particles interacting via the potential specified in Eqn l2.H the phase 
diagram is not known even when the external laser field is absent . However, Brownian 
dynamics simulations for a similar system 19] suggest that at Uq ~ 3.1, there is a transition 
from a single, dispersed phase to a phase where the colloidal particles start forming clusters. 
In our present work, we consider the case Uq = 5.4 and r\ = 30. The first parameter defines 
the depth of the attractive potential and the second one fixes its width (cx/30). In the 
range of densities of our interest, 0.85 < p < 0.95, Bolhuis et al j20j have observed the co- 
existence of a high-density solid with a dilute gas for a two-dimensional system of particles 
interacting via a short-range attractive square potential. Co-existing gas and solid phases 
are also found in three dimensions |2lJ when the inter-particle attraction has a short range. 
Our simulations also show the formation of clusters for the choice of Uq and rj mentioned 
above. In FigJTJ we have plotted the average density, representative of our simulated system, 
at p = 0.86. The plot clearly shows the co-existence of crystalline clusters with a gaseous 
phase, with particles arranged in a triangular lattice inside the clusters, as observed in the 
Brownian dynamics simulations jl^. The structure factor, S(q), for a cluster, plotted in 
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Fig |21 shows peaks corresponding to a triangular lattice structure with spacing equal to a, 
the hard sphere diameter. 

Our objective is to study how such a system of colloidal particles behaves in the presence of 
an external potential with one-dimensional periodic modulation d which is not commensurate 
with the distance between the lattice planes of the triangular structure that the particles 
form in the absence of the field. For this purpose we calculate the average density p(x, y) and 
the structure factor S(q) for different strengths of the potential at a fixed particle density. 

For p = 0.86, when the strength of the potential is low, our simulations show that the 
system continues to form crystalline clusters with inter-particle separation a. In Fig. Efa) 
and FigEJ^a) , we have plotted respectively the average density and the corresponding S(q) 
for the colloidal system when the strength of the laser-induced potential is f3V e = 0.50. As 
can be seen from the plots, the nature of S(q) is similar to the case where the potential is 
absent. 

When the strength of the potential is increased to (3V e = 2.0, the clusters break up and the 
particles try to align themselves along the potential troughs. The average density p(x, y) 
for this value of (3V e , plotted in FigJH^b), shows that although there is a local triangular 
structure, the particles are also getting arranged in the y direction, due to the influence 
of the external potential. The 5(q) for this situation, plotted in FigHfb), shows that the 
height of the peaks corresponding to the triangular lattice has decreased considerably com- 
pared to that in FigHf a) and new peaks have begun to appear on the q y axis with q y values 
corresponding to the wave- vector of the applied field. This clearly is an effect of the incom- 
mensuration. The system is in a 'frustrated' state - the particles are making an attempt 
to lie at the troughs of the periodic potential, but the energy they gain by doing so is not 
sufficient to break the 'bonds' of the triangular lattice with spacing a. 

In Fi^Sfc), we have plotted p(x, y) for the system when the strength of the laser field has 
been increased to (3V e = 5.0. The crystalline clusters have now melted - the particles have 
become confined along the potential troughs. The S(q) for the system, plotted in Fi^fc), 
is characteristic of a modulated liquid, with the peaks of the structure factor located only 
at multiples of the characteristic wave- vector q y = 2ir/d of the external periodic potential. 

When the strength of the periodic potential is increased to higher values, the motion 
of the particles in the direction transverse to the potential troughs decreases considerably. 
The average density for the particles at a potential strength of f3V e = 100.0, plotted in 
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FigEtd), shows crystalline order corresponding to a square lattice, which is also reflected 
in the corresponding S(q), shown in FigEfd). The peaks of the structure factor on the q y 
axis occur at wave vectors that correspond to multiples of the wave-vector of the external 
potential, whereas the peaks on the q x axis correspond to an inter-particle separation of a 
inside each potential trough. 

We have also computed the spatial correlation function g(x), which is the pair correlation 
function along the potential troughs, for the two phases observed for the potential strengths 
/3V e = 5.0 and (5V e = 100.0. The maxima of the correlation function for (3V e = 5.0, plotted in 
Fi^Hfa), decay quite fast, and at long distances g{x)— >1. This is characteristic of liquid-like 
behaviour for the particles confined in the potential troughs. However, when j3V e = 100.0, 
from the g(x) data plotted in Fij^b), we can conclude that there is a periodic modulation 
of the density of the colloidal particle in the potential troughs, the period of this modulation 
being a as expected. 

To check whether the formation of the crystalline solid at the higher values of (3V e is 
a finite size effect, we have calculated the corresponding order parameter p q for different 
system sizes, N = 1600, 1024, 900. For the solid phase, the order parameters are the Fourier 
components of the density p(r) calculated at the reciprocal lattice points {q}. From the 
peaks of the structure factor plotted in Fig0Jd), we can get the four smallest reciprocal 
lattice vectors for the square lattice formed by the particles for f3V e = 100.0. Of these four 
vectors, the two lying on the q y -axis correspond to the ordering due to the external field 
and the p q for these vectors are also non-zero for the modulated liquid. So the p q for the 
other two wave- vectors, denoted by p q , are the relevant order-parameters for checking the 
effects of the finiteness of the system. The average order parameter < p^ >, has the values 
0.731, 0.663, 0.661 for iV = 900, 1024, 1600 respectively, implying weak dependence on the 
system size. Therefore, we can conclude that the square lattice formed at [3V e = 100.0 
represents a crystalline phase. 

The appearance of this crystalline phase at a large value of (3V e can be understood from 
the fact that, for p = 0.86, the inter-trough spacing, which is governed by the periodicity 
of the external field d, is approximately equal to a, the distance at which U(rij) has its 
minimum. In fact, for 0.811 < p < 0.866, the spacing d is within the range of the attractive 
part of the potential, 1.033 < r/a < 1.0, for the value of rj used in our simulations. Thus, 
at p = 0.86, if the transverse motion within a trough is suppressed, correlations develop 
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across the troughs resulting in ordering of the particles in the y direction, in addition to 
the ordering with spacing a along the troughs. Such a crystal structure corresponds to the 
lowest energy configuration consistent with the density and the inter-trough spacing. 

Thus the system of colloidal particles, which had preferred to form crystalline clusters at 
low field strengths, regains a different crystalline phase at high field strengths, after passing 
through an intermediate modulated liquid phase. This phenomenon may be called re-entrant 
freezing. It is interesting to note that the two crystalline phases have different symmetries, 
which is a consequence of the incommensuration effects. Another aspect to note is that 
the solid is in co-existence with a dilute gas and the 'voids' corresponding to the gaseous 
phase can occur at any point in space. This is a probable reason for the decrease in the 
height of the maxima of g(x) at large distances. The inter-trough correlations, observed at 
(3V e = 100.0, are not developed at lower field strengths because the motion of the particles 
in the transverse direction is not suppressed sufficiently and the particles can vibrate within 
the width provided by the trough, thereby causing a destruction of the crystalline order. 

However, we should note that the kinetics of systems with hard-core interactions may 
depend strongly on the initial state j^jj. In our simulations, we observe such behaviour. For 
some initial configurations, the system gets kinetically jammed and due to this jamming, 
even at high field strengths, there is a co-existence of the square and triangular phases. In 
simulations, one can escape from jamming if a combination of local and non-local Monte 
Carlo moves are used. 

For p > 0.866, the wavelength of the external potential becomes smaller than the particle 
diameter and therefore, for these densities, the square lattice become unfeasible at high 
potential strengths. However, even for these densities, the particles in neighbouring troughs 
would like to be in contact with each other and that can only happen if they form a rhombic 
structure. In FigEJ we have plotted the average density for p = 0.92 at j3V e = 100.0. The 
particles have formed a 'mixed crystal' , i.e they have formed domains which have local 
crystalline structure similar to the two structures shown in the figures below the p(x, y) 
plot. As (3V e — > oo, the motion of the particles becomes one-dimensional and the system 
tries to attain one of these two crystalline structures in order to minimize its free energy. 
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B. Model 2 



For the other choice of the pair potential given in Eqn l2.4l we have considered 400 particles 
in a rectangular box commensurate with a triangular lattice structure and with periodic 
boundary conditions. The screening length is fixed at na s = 15.0 where a s , the lattice 
spacing in a triangular arrangement, has been used as the unit of length in all the expressions. 
Without any external field, the system (colloidal particles of diameter 1.07/im, surface charge 
Z* = 7800e and density n p = 1.81 x 10 7 / cm 2 , suspended in water having dielectric constant 
e = 78 at a room temperature of 298K) is in the liquid phase. The parameters for the 
Gaussian well are chosen as A = ^ and B = 50. In the first part of our simulations for this 
system, the position of the attractive minimum is set at A = 1.3. 

To observe the effect of the laser field in this system of particles, the average density is 
calculated for different values of j3V e . In Fig[?l we have shown a contour plot of the average 
density, representing a section of the simulation box, in the absence of the external potential. 
The average density does not show any order, signifying a liquid state. 

FigJS] contrasts the effect the external potential has on the average density depending on 
whether the inter-particle potential has the attractive part or not. At j3V e = 0.4 (subplots 
(a) and (b)), the liquid-like structure of Fig|7|is replaced by one corresponding to a triangular 
lattice structure in both cases. The order parameter p^, corresponding to the triangular 
lattice structure, has the value 0.6 for the pure DLVO potential and 0.53 for the case when 
the attractive part is present. As the field is increased to /3V e = 1.5, the crystalline order for 
the pure DLVO potential becomes much sharper, see Fig|Sfc), and the value of p^ (for the 
triangular lattice) becomes higher. But with the attractive part present (see FigJHfd)), the 
average density looks like a set of liquid-like strings, indicating a modulated liquid phase, 
with Pq(for the triangular lattice) becoming much smaller « 0.15. At still higher fields, 
/3V e = 4.0, nothing changes qualitatively for the pure DLVO case (see Fig|Hfe)). However, 
the average density plot for the case with the attractive potential (FigJHff )) indicates the 
formation of a rectangular lattice. 

Therefore, for this choice of inter-particle potential, we again find that the system of 
colloidal particles, which had formed a crystalline structure for low strength of the exter- 
nal potential, melts as the value of j3V e is increased. But at stronger external potential, 
the particles again form a crystalline structure which has rectangular symmetry, similar to 
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what we had observed for the case of Model-1. Hence, here also, we observe a re-entrant 
crystallization. 

At high external potential strengths, different choices of the parameter A need not always 
result in a rectangular lattice. In Fig. E]we have plotted the average density for various other 
values of A and f3V e = 4.0. It shows that a rich variety of phases can arise. To understand 
the structures, we consider the potential energy of a given particle in either triangular or 
rectangular lattice arrangement. Since both the DLVO and the Gaussian attractive part 
of the potential fall off rapidly with increasing particle separation, the contribution from 
the first shell of neighbours is the most significant. Because of the large f3V e , the potential 
rises sharply in the y direction from the potential troughs. Thus we look for the potential 
experienced by a particle along the x direction with the y co-ordinate fixed to be at a 
minimum of the external potential. Considering the six nearest neighbours for the triangular 
lattice case and the eight nearest and next-nearest neighbours for the rectangular lattice, 
we have calculated the effective potential well when all the neighbours are in perfect lattice 
positions (Fig. HT^i . For A = 0.8, the average density (Fig. Eta)) looks like a modulated 
liquid with some superposed modulation perpendicular to the field direction also. Fig. HOf a) 
shows the behaviour of the potential well for the same parameter values - both triangular 
and rectangular lattices are energetically unstable. The energy minimum is for a hexagonal 
arrangement of the neighbours similar to that in triangular lattice, but the central particle 
has two equivalent energy minima displaced from the center. Thus the system remains 
frustrated with some residual hexagonal order. 

As A is tuned to a value of 1.0, a deep minimum in the single particle potential cor- 
responding to a triangular lattice develops (Fig. ITUT b^. However, since we increased A 
through a frustrated potential structure, the system takes a long time to relax to an unique 
triangular lattice spanning the whole system. Fig. |^c) shows the average density plot for 
A = 1.0. To accommodate the history of the frustrated structure at lower A, the shown 
part of the system goes through a coherent shift along the x direction. At long times, the 
average density shows sharp contours corresponding to a triangular lattice. 

With increasing A, the single particle potential with the neighbours in a regular lattice 
structure again shows degenerate minima, but this time the rectangular lattice wins over 
the triangular lattice energetically. For A = 1.2 the single particle potential is shown in 
Fig. ITUf c) and the corresponding average density contours in fig. Ete). The structure is that 
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of modulated liquid, with superimposed rectangular modulation. 

When A = 1.3, the rectangular structure is energetically stable and more favoured than 
the triangular lattice (Fig ITUT d)). The average density plot (Fig. [HIT)) shows sharp contours 
corresponding to a rectangular structure. 

Therefore, once again, we observe that as a function of the location of the minimum 
of the attractive part of the inter-particle potential, the system switches from one kind of 
crystalline structure (triangular) to another ordered structure (rectangular) via a modulated 
liquid phase. 

IV. CONCLUSIONS 

In summary, in this paper we have investigated the effect of an external laser field, 
periodically modulated in one dimension, on a system of colloidal particles confined to two- 
dimensions and interacting via a potential that includes a short-range attraction component. 
The presence of this attractive interaction introduces a new length scale (corresponding to 
the minimum of the attractive pair potential) which can be incommensurate with the wave- 
length of the potential due to the externally applied laser field. We find that the competition 
between the two incommensurate length-scales results in the phenomenon of reentrant crys- 
tallization at high field strengths. In both the models we have studied, we find that the 
crystalline phase attained by the particles at low field strengths melts into a modulated 
liquid when the strength of the field is increased. However, it eventually crystallizes again, 
at higher field strengths, due to suppression of the particle motion in the direction of mod- 
ulation of the applied potential. This phenomenon is opposite to the situation when the 
attractive interaction is absent, where one observes only re-entrant melting. We have also 
observed, at least in one of the models, that depending upon the position of the minimum 
of the attractive potential, one can get different phases when both the particle density and 
the laser potential are kept constant. 

We hope that our simulations will motivate experiments to observe the effect of the laser 
field in two-dimensional colloidal mixtures. Also, in the context of colloidal aggregation, 
our simulations show that structures with square symmetry in the aggregated phase can be 
stabilized in the presence of substrates formed by laser intensity fringes. We have also shown 
that other novel structures can be formed by tuning various parameters of the colloidal sys- 
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tern. Further theoretical studies to understand the phenomena observed in our simulations 
in greater detail would be interesting to pursue. 
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FIG. 1: Model 1. Average particle density p(x,y) at p = 0.86, for the case when the external field 
is absent, showing co-existence of crystalline clusters with a gaseous phase. 
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FIG. 2: Model 1. Struture factor 5(q) of a colloidal cluster at p = 0.86, for the case when the 
external field is absent, showing the existence of triangular lattice structure within the cluster. 
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FIG. 3: Model 1. At p = 0.86, the average particle density p(x,y) shows (a) 2d triangular lattice 
when f3V e = 0.5, (b) 'frustrated' liquid when (3V e = 2.0, (c) modulated liquid when (3V e = 5.0, (d) 
2d square lattice (3V e = 100.0. Thus, at this density, the system undergoes re-entrant freezing as a 
function of the strength of the externally applied field. 
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FIG. 4: Model 1. For p = 0.86, Struture factor S(q) shows (a) 2d triangular lattice when (3V e = 0.5, 
(b) 'frustrated' liquid when (3V e = 2.0, (c) modulated liquid when j3V e = 5.0, (d) 2d square lattice 
when pv e = 100.0, confirming the re-entrant freezing experienced by the system as (3V e is increased. 
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FIG. 5: Model 1. Pair correlation functions g(x), at p = 0.86, along the potential troughs for (a) 
modulated liquid (when (3V e = 5) and (b) square lattice (when (3V e = 100). 
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FIG. 6: Model 1. (a). The average density for p = 0.92 at (3V e = 100, showing that the particles 
form a mixed crystal. 

(b). The two panels are snapshots of the crystalline structures having same energy, one of which 
will be formed by the particles as (3V e — > oo, at p = 0.92. 
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FIG. 8: Model 2. The left hand panel (a, c and e) shows the average density for simulation with 
just the DLVO potential. The right hand panel (b, d and f) shows the same for simulation with 
attractive potential. 
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FIG. 9: Model 2. Average density for na s = 15, (3V e = 4 and with different values of A. 
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FIG. 10: Model 2. Potential experienced by a given particle with neighbours in perfect hexagonal 
(points) or rectangular (lines) arrangement for different values of A. 
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